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Abstract Visualization using tracer particles is a relatively new tool available for the study of 
superfluid turbulence and flow, which is applied here to oscillating objects submerged in the liquid. We 
report observations of a structure seen in videos taken from outside a cryostat filled with superfluid 
helium at 2 K, which is possibly a vortex loop attached to an oscillator. The feature, which has 
the shape of an incomplete arch, is visualized due to the presence of solid H 2 tracer particles and 
is attached to a beam oscillating at 38 Hz in the liquid. It has been recorded in videos taken at 
240 frames per second (FPS), fast enough to take ^ 6 images per period. This makes it possible 
to follow the structure, and to see that is not rigid. It moves with respect to the oscillator, and its 
displacement is in phase with the velocity of the moving beam. Analyzing the motion, we come to 
the conclusion that we may be observing a superfluid vortex attached to the beam and decorated by 
the hydrogen particles. An alternative model, considering a solid hydrogen filament, has also been 
analyzed, but the observed phase between the movement of the beam and the filamentary structure 
is better explained by the superfluid vortex hypothesis. 

Keywords Quantum fluids • Superfluid Helium • Flow visualization • Vortex loops 
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1 Introduction 


In the superfluid phase of liquid Helium, below 2.177 K, the circulation is quantized in units of a 
flux quantum (k = h/rriHe — 10 ~ 3 cm 2 /s where m# e is the mass of a He 4 atom and h is the Planck 
constant). The existence of vortices with a single flux quantum was independently proposed by 
Feynman and Onsager 1|2] and the first measurements showing quantized circulation were made by 
Vinen (3]. More recently superfluid vortices have been observed by Bewley et al [4] , and this group has 
developed solid hydrogen tracers [5|6] to visualize the flow and has observed many interesting features 
of vortex physics such as re-connections [4j[5] and Kelvin waves [T] . Visualizations of turbulence 
generated by counterflow have also been obtained by this technique [8] 10 which is becoming a 
powerful tool in the study of Quantum Turbulence 111 . 

However, apart from some preliminary work 12 , the visualization of flow around objects oscillat¬ 
ing in superfluids has not been explored so much, although Vinen and Skrbek [13,14 have pointed 
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Figure 1 a)Diagram of the vibrating beam, with measurements in millimeters. The position of the driving magnet 
is shown as a circle in this frontal view of the setup, b) Position of the points x and xo (see text) in the approximate 
perspective seen in the videos. Vq represents the local superfluid velocity due to the heat load of the cryostat. 
The driving magnet and its position are shown as a cylinder in perspective, c) Proposed geometry of the vortex 
loop and definition of some quantities used in the analysis d) Alternative mass-spring model considering a possible 
hydrogen filament, instead of the loop shown in c). 


out that tracer imaging of superfluid oscillatory flows with a classical analogue could provide valuable 
information. We have recently developed a simple system using solid hydrogen particles for visual¬ 
izing flow around oscillating objects in superfluid helium 12 . In the following we present further 


observations made using this system which show a behavior which is consistent with the presence of 
a superfluid vortex half-loop attached to a beam oscillating in liquid helium. The loop is observed 
to expand and contract and is attached to the beam throughout our observation. To our knowledge, 
such behavior has not been directly visualized previously, although vortex loops have been observed 
indirectly through their attachment and detachment 1 15 j and discussed theoretically before 16-18 . 
We therefore believe that our observations could shed new light on the problem of superfluid tur¬ 
bulence, in particular on the behavior of vortices attached to a solid boundary with oscillating flow 
and their stability and pinning. 


2 Experimental Details 

The experimental arrangement has been described in detail elsewhere [12 . The system studied 
consists of a vibrating beam, driven magnetically by using a permanent magnet attached to the 
beam and a coil fixed to a rigid frame. In this way the beam oscillates with velocity perpendicular to 
its wide dimension at a frequency of 38 Hz. A sketch of the geometry is found in Fig. [l] Videos are 
taken with a camera at 240 frames per second (FPS) so the time interval is 4.17 ms between frames 
and we use this as our time reference to calculate velocities. The helium temperature was 2.07 K 
throughout the experiment. 

A mixture of one part hydrogen to 50 parts helium gas at 500 torr is introduced from room 
temperature to form the solid hydrogen particles and around a hundred cubic centimeters of gas are 
injected each time. To illuminate the tracer particles we used a green laser beam. The frozen H 2 
particles are not expected to absorb significant energy in the visible [19]. The laser is on the outside 
of the dewar and the light passes through an optical fiber which ends less than a centimeter away 
from the oscillating beam, illuminating the particles perpendicular to the line of sight of the camera. 
The fiber is polished at the end, giving a three dimensional cone of light. Distances in the image are 
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calibrated with respect to the measured dimensions of the small magnet (a cylinder 5mm long and 3 
mm in diameter), which is used to drive the oscillator. The size of a pixel corresponds to around 70 
microns in the object but the light could be scattered from particles which are smaller than this. It is 
hard to evaluate the minimum observable dimension, but we estimate our particles to be distributed 
in size from well below 70 microns to 200 microns. 

An important feature of our setup is that we are forced to remove the outer nitrogen Dewar to 
avoid the blurring of the images produced by the boiling nitrogen. For this reason, the heat load is 
considerable. We can estimate it by measuring the volume of He evaporated as a function of time 
and using the known latent heat of evaporation. The heat input Q is not constant, but has been 
roughly measured to he between 800 and 470 mW. The helium Dewar is 6 cm in diameter, so the 
calculated counterflow velocity vq is between 0.085 and 0.05 cm/sec if we assume a uniform heat flow. 
These numbers are only rough estimates since the geometry is not simple, the heat input comes from 
radiation through the walls, conduction down the glass walls, etc. On the Dewar there is a flange of 5 
cm diameter about 5 cm above the vibrating beam which also complicates the counterflow geometry. 


3 Experimental Results 

The main observation of the experiment is the existence of a structure formed by the H 2 particles 
that seems to be attached to the beam and oscillates with it, though not in a rigid fashion. It seems 
to elongate and contract when the beam oscillates, and has a curved shape, somewhat resembling 
an incomplete arch. This structure is seen in all videos taken during the experimental run. We show 
still images in Fig. [2] taken over one complete cycle of the oscillator. The images are amplified close 
to the maximum resolution and the pixel structure can be clearly seen. 

For analysis, we have chosen to follow two points marked in Fig. [2] and Fig. [l] as x and xq. Point 
Xq is fixed to the oscillating beam and x is a point on the arch, which is seen to shrink and grow 
periodically. We have followed x and Xq over several cycles of the oscillator, observing the images by 
eye, and digitizing the positions x and xo by means of a computer. The difference in the position of 
these two points, which have a sinusoidal motion, is proportional to the length of the arch formed 
by the decorating particles. We have taken a well defined and easy to follow point for xq instead of 
the base of the arch, so x — x$ has a constant displacement superposed to the periodic component 
but the sinusoidal variation is proportional to the length of the arch. 

Digitizing the positions of x and Xq we are able to calculate the velocities of the beam and the 
arch as a function of time (Using the x$ positions in successive images and dividing by the known 
time of 4.17 ms between frames). We can also evaluate the distance (x — Xq) in each frame. The 
results are shown on Fig. [3] Open circles correspond to the velocity of the beam and filled circles 
to the distance x — xq. A few cycles are shown, and they include information from three different 
video sections, all taken the same day but at different times. We have fitted the experimental points 
with sinusoidal functions by least squares, and the results are shown as lines in the plots. The fits 
show clearly that the velocity of the beam is in phase with the relative displacement x — xq. Bearing 
in mind that x — Xq is proportional to the length of the arch formed by the solid particles, the 
length of this structure therefore is in phase with the velocity of the beam. We have also evaluated 
the velocities of x and Xo separately, and find that they are around 90 degrees out of phase with 
each other. Therefore the arch structure is not rigidly attached to the beam, but shows an internal 
shrinking and growing movement in phase with the velocity. 

The phase difference between velocities and x — Xq is shown in a different way in Fig. [4j We show 
the data in a parametric plot, on the left hand side the velocity of the the beam ( Vseam ) is plotted 
against a point corresponding to the top of the arch (Vx) and it is seen that the plot is almost 
circular, as corresponds to a parametric representation of two sinusoidal functions with a 90 degree 
phase difference, while on the right a plot of x — Xq against VBeam shows a very elongated ellipse at 
a 45 degree angle, as would be seen for two sinusoidal functions that are almost in phase. 

Although we only show a few periods, videos are longer than this, and the motion of the arch was 
observed to remain basically unchanged throughout the whole experimental run. The video camera 
was not always on, but the oscillation of the beam was not changed for 13 minutes. During this time 
we filmed ten sections of video at 240 and 480 frames per second, and in all of them we see the loop, 
with similar behavior. The sections shown are representative of the first and last parts filmed, at 
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Figure 2 Still images of a video. The six images correspond to successive frames, taken at 4.17 ms intervals and 
therefore covering a complete cycle of the oscillator whose frequency cjb/ 27t = 38Hz. The pixel size can be seen at 
this magnification. 


240 FPS. We have not included data for the 480 FPS because the resolution and lighting are not 
good enough for measuring quantitatively, although the images show movement of the loop that is 
compatible with that seen at 240 FPS. Several experiments with the same setup were performed, 
but the structure seen here was seen in only in this particular run when the videos were analyzed in 
detail later. It seems therefore that the formation of a loop is not a reproducible feature, but depends 
on several uncontrolled factors as is expected in turbulent regimes. 


4 Discussion 


We have come to the conclusion that the behavior observed is consistent with the motion expected 
from a vortex half-loop pinned to the beam. This appears to cover the main facts, although other 
possibilities have also been considered. We use for our analysis the relationship obtained by Schwarz 
[201 for a stable vortex loop. According to his (20 Eq. 20 a vortex loop of radius Ro will not shrink 
or grow if it moves with a velocity v s with respect to the superfluid: 


v s 


hi 8 Rq 

AnRo _e 1 / 4 ao 


(i) 


with ao ~ 1.3 x 10 -8 cm an adjustable parameter roughly equivalent to the size of the vortex core. 
Conversely, if the loop is fixed in position, a flow of the superfluid with velocity v s maintains a stable 
radius, if the sign of the velocity and the vorticity of the loop are in the right orientation. Section 
IV A of Ref. (20| is also relevant to our situation, since it discusses a bent vortex loop pinned at two 
points and his Fig. 29 shows the shape calculated for different values of superfluid velocity. In our 
case, using Eq. [l]and an average curvature radius of 0.3 mm the velocity v s for a stable half loop 
would be 0.027 cm/s. The estimate for the average counterflow velocity Vq is a factor between 2 to 3 
times greater than v s . However the flow due to Vq close to the beam would probably be lower than 
Vq itself due to the obstacles present. If we assume that Vq along the beam is modulated by the 
periodic motion of the beam, we could explain the stretching and contraction of the arch attached 
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Figure 3 Comparisons of beam velocity Vseam and x — xo for three sections of the videos. Open circles correspond 
to the velocity of the beam and filled circles to the distance x — xo. The dashed and full lines are sinusoidal fits to 
the data and it can be seen from the graphs that Vseam and x — xq are almost in phase. 


to the oscillating beam. The modulation due to the beam Vseam and the counterflow velocity would 
add to produce a time dependent velocity 


VOsc(t ) = A • Vq + B ■ V B eam Sm(ut) 


( 2 ) 


where A and B are adjustable parameters, Vseam is the measured velocity of the beam, and uo the 
oscillating frequency. 

From Eq. [l] we can obtain an approximate expression for the radius of the loop due to the time 
dependent velocity 

4:7t(A Vq B V B eam SU1 (ut)) ^ ^ 

Where Rq = A7T ^. Vq * C ~ 0.3 mm would be the stable radius if Vq has no modulation, and we 
have not considered significant the modulation of the logarithmic term, including it as a constant 


C = In 


SRp 


2 1 / 4 Q 


We have fitted (x — xq) ~ R(t ) to the expression given by Eq. 


and we show the 
37 and B = 0.001, 


results in the bottom panel of Fig. [4j The adjustable parameters used are A — C 
and a shift has been introduced to compensate for the arbitrary origin when choosing xq. In this 
case we have used the lower value of Vq = 0.05 cm/s. The fit is quite good, and although it is not 
a least squares fit and the parameters were chosen by hand, the equation is capable of reproducing 
the observations. 

A second possibility is that the structure seen is not a vortex loop, but a filament of solid hydrogen, 
as has been observed by Gordon et al |2lJ 22 . We could model a solid filament as a mass attached 
to a spring, as is shown in Fig. 0d). The movement of the attaching point Xo would move the mass 
at x as a forced harmonic oscillator 


mx + 7 ^ + k(x — Xq) = 0 (4) 

the well known solution for this equation is harmonic motion, with a phase difference 0 between 
and x. Our observed phase difference is of around 90 degrees, between these quantities, which would 
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Figure 4 Upper panel: Parametric plots of Vx vs Vseam (left), and x — xo vs Vseam (right). Open circles: 
experimental points, lines: sinusoidal fits to the data, equivalent to the full and dotted lines shown in Fig. 3. The 
almost circular plot on the left indicates that the corresponding sinusoidal motions are 90 degrees out of phase, 
while the highly tilted ellipse at the right implies that the oscillations are almost in phase. Lower panel: Fit of 
x — xo oc R(t) using Eq. [3] and fitting parameters A = 0.37 , B = 0.001. Lines and open circles correspond to the 
measured x — xo and the fit is shown in dotted lines. The good agreement indicates that the vortex loop model 
gives an adequate description of the observations. 


be the case only if the driving frequency ujb/ 27t of 38 Hz were accidentally close to the resonance 
frequency of the filament ujf = yjk/rn. It is highly unlikely that ujf is the same as the frequency of 
the beam ujb, and it could be expected that </> would either be zero (if cj/ » uob) or 180 degrees (if 
ujf « ojb)- A third possibility is that the oscillating filament could be overdamped by the influence 
of the normal component, but in this case x would follow the fluid, which in the reference frame of 
the laboratory, moves 180 degrees out of fase with the point Xo belonging to the beam. In fact, our 
observation could be a combination of a hydrogen filament and part of a vortex loop, attached to 
the end of the filament and a point in the beam. The filament shown in Fig. 3 of Ref. 22 shows 
movement in the normal fluid, but to have the phase observed here some form of vortex section, 
closing the loop and with movement given by Eq. [3] seems to be necessary to explain the observed 
behavior. 

For further analysis, the local approximation could also be used, and in this approach the dy¬ 
namics of a quantized vortex is described by the equation proposed by Schwarz 23 


ds 

— = /3s' x s" + v s + as' x (v n - v s - /3s ' x s") 


( 5 ) 


here s is a point on the core of a vortex loop, /3, a the coefficient of mutual friction, v n and v s the 
velocities of the normal and superfluid fractions, s' is the tangent to the vortex core, s" the principal 
radius. The proposed structure of the vortex half-loop and the definition of the vector quantities s' 
and s" in Eq. [5] are shown in Fig. [lj We can assume that the vibrating beam pushes both the normal 
fluid and the superfluid together, as well as modulating Vq as described above, so we have 


v n • r = v s • r = v\ sin ujt • cos p (6) 

with p the angle between the local velocity and the vortex. Then Eq. [5] implies 

dR a/3 

—— = D • v\ sin c ut • cos cp -— 


( 7 ) 
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The equation is local, so that ^ changes over the circumference with cos p. It also changes the 
shape of the half loop, depending on the relative orientation of loop and velocity, but we can get an 
approximate value for the average radius (R) neglecting the second term with respect to the first 
and integrating in time and over p 

( R ) = Rq — D— cos uot. (8) 

uj 

with D a parameter taking into account the angular integration over p. We do not have enough 
resolution to detect the changes of shape implied by Eq. [8j but it would appear that the effect is 
smaller than that of Eq. [3j Furthermore it gives motion with a phase that is at 90 degrees from the 
velocity of the beam, instead of the in phase motion observed. 

In fact, Eq. [3] can be taken as a quasi static non local solution including a modulation of Vq, and 
Eq. [8] as a local time dependent correction due to the presence of v\ whose importance is given by 
the parameter D. The good fit obtained with Eq. [3j seen in Fig. [2] seems to indicate that the effect 
of the correction of Eq. [8] is not very large, although it could be responsible for the fact that in the 
measurements x — Xq and Vseam are not always in phase. In fact, seen over many cycles, there are 
small irregularities in the motion, where the structure seems to stretch more or less. However, the 
overall stability is preserved, as stated earlier, over at least the 13 minutes where we have partial 
observations. Since the frequency is 38 Hz, the number of cycles over which the stretching and 
shrinking is repeated is of order 3 x 10 4 . 

In conclusion, we have observed, by decoration with solid H 2 tracers, a structure which moves 
attached to a vibrating beam. From an analysis of possible models for the observed motion, we 
conclude that it behaves as expected for a vortex half-loop attached to the oscillator. This accounts 
for the phase relationship between position and velocity, and we obtain a good fit between the model 
and the video images, while an alternative explanation postulating a hydrogen filament requires 
an unlikely coincidence between the driving frequency and the natural frequency of the hydrogen 
filament. A third model, considering that the motion is due to the drag of the normal component on 
an over damped filament, would also produce a phase difference that is not the one observed. 
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